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The  AURORA  simulator  is  designed  to  provide  short  pulses  of 
high  energy  bremsstrahlung  over  volumes  large  enough  to  irradiate 
entire  weapon  systems  or  large  subsystems  uniformly.  Typically, 
120-ns  pulses  of  20  to  50  krads(Si)  can  be  delivered  into  a 1-m 
cube.  Because  of  the  large  irradiation  volume,  it  is  not  unusual 
for  an  experimenter  to  desire  dose  measurements  at  up  to  200  loca- 
tions within  his  system,  sometimes  with  a pulse  repetition  rate 
of  four  pulses  per  day.  Most  frequently,  the  dose  data  are 
requested  in  rads (Si). 

A dosimetry  system  for  this  facility  therefore  required  the 
following  characteristics: 

(a)  Ability  to  measure  equilibrium  and  nonequilibrium 
silicon  dose 

(b)  Useful  dynamic  range  extending  from  1 to  10^  rads (Si) 

(c)  Independence  of  dose  rate  up  to  5*10^^  rads (Si) /s 

(d)  Capability  of  reading,  processing,  and  recording  up 
to  800  data  points  per  day. 

This  report  describes  a dosimetry  system  designed  to  fit 
these  requirements.  The  rationale  for  our  choice  of  thermolumi- 
nescent dosimeters,  characteristics  of  these  dosimeters,  design 
of  instrumentation  and  data  analysis,  precision  and  accuracy 
are  described. 
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1 . INTRODUCTION 


AURORA,  a generator  of  intense  bremsstrahlung  pulses,  is  used 
principally  for  simulating  the  effects  of  a nuclear  weapon  prompt  gamma 
pulse  on  military  electronic  systems.  AURORA'S  unique  feature  is  its 
extremely  high  volumetric  dose  output;  it  can  deliver  25  krads(Si) 
throughout  a 1 m-diam,  1 m-long  cylindrical  volume,  or  50  krads(Si) 
throughout  a 25-cm  sphere.  AURORA  is  located  at  Adelphi,  MD,  and  is 
operated  by  the  Harry  Diamond  Laboratories  for  the  Defense  Nuclear 
Agency. 

This  report  describes  the  techniques  and  instruments  used  at  AURORA 
for  measuring  the  radiation  dose  and  dose  rate.  This  is,  of  course, 
important  to  experimenters  at  AURORA  for  understanding  and  interpreting 
their  dose  data.  But,  more  than  this,  it  is  hoped  that  this  information 
will  be  useful  also  to  anyone  who  must  do  dosimetry  at  the  many  other 
pulsed  bremsstrahlung  facilities  that  now  exist  or  are  planned. 

The  present  report  describes  in  detail  the  experimental  techniques 
of  dosimetry  as  they  are  now  implemented  at  AURORA.  It  thus 
complements,  but  does  not  significantly  reiterate,  the  contents  of  two 
other  reports  concerned  with  AURORA  dosimetry.  The  first ^ describes  the 
experiments  and  rationale  that  led  to  the  eventual  adoption  of  the 
AURORA  dosimetry  system  .as  it  now  exists.  It  provides  especially 
valuable  insight  into  the  characteristics  of  the  various 
thermoluminescent  phosphors  and  their  applicability  to  the  AURORA 
dosimetry  problem.  Some  of  this  information  can  be  found  also 
elsewhere.^  Another  report  (unpublished)  does  not  dwell  on  technique 
per  se,  but  describes  the  radiation  output  characteristics  of  AURORA  to 
the  extent  that  they  are  known  at  this  time.  These  three  reports 
together  tell  the  whole  AURORA  dosimetry  story. 

Readers  who  are  interested  in  the  general  characteristics  of  the 
AURORA  simulator  and  facility  operations  and  instrumentation  should 
consult  the  AURORA  Facility  brochure  (available  from  author) . 

2.  DOSIMETRY  THEORY 

This  chapter  is  not  intended  to  be  a definitive  treatise  of 
dosimetry  theory.  The  aim  is  to  present  just  enough  theoretical 
background  here  to  justify  some  of  the  concepts  and  equations  used  in 
the  following  chapters.  Readers  desiring  a more  rigorous  treatment  of 

^F.  H.  Attix,  Dosimetry  for  the  Aurora  Facility , Final  Report,  Naval 
Research  Laboratory,  Washington,  DC,  NRL  Memorandum  Report  2552,  Revised 
(March  1973) . 

^S.  G.  Corbies,  F.  H.  Attix , and  K.  Kerris,  Thermoluminescent 
Dosimeters  for  High-Dose  Applications,  Health  Physics,  2^  (1973),  499. 
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dosimetry  theory  are  urged  to  consult  other  sources  for  general 
background  and  a discussion  of  charge  particle  ecjuilibrium^ ' ^ and  for 
Bragg-Gray  cavity  theory. 

Radiation  dosimetry  is  the  measurement  of  the  energy  deposited  in 
matter  by  radiation.  The  energy  per  unit  mass  deposited  in  an 
infinitesimal  volume  of  a material  is  called  the  dose.  That  is,  if  we 
let  be  the  energy  deposited  in  a very  small  volume  of  a substance  of 
mass  Am,  the  dose  is  given  by 


D 


Am 


(1) 


The  unit  of  dose  is  the  rad,  which  is  equal  to  100  ergs/gram  by 
definition.  Since  the  quantity  of  energy  absorbed  from  any  given 
radiation  field  depends  on  the  nature  of  the  material,  it  is  necessary 
to  specify  which  particular  material  one  is  referring  to  when  talking 
about  dose.  This  reference  is  always  made  by  a parenthetical  notation 
immediately  following  the  name  of  the  units?  e.g.,  rad(air).  Thus,  if 
exposure  to  a certain  source  of  radiation  results  in  600  ergs  being 
absorbed  by  1 gram  of  silicon,  we  say  that  the  dose  was  6 rads (Si). 


In  the  particular  situation  of  interest  at  AURORA,  that  is,  where 
the  radiation  consists  of  energetic  photons,  energy  is  deposited  in 
matter  in  a two-step  process-  First,  a photon  transfers  energy  to  an 
electron  in  the  material  by  the  photoelectric,  Compton’,  or 
pair-production  process-  Then  this  electron  transfers  its  kinetic 
energy  to  the  material  via  ionization  processes. 


The  dose  deposited  in  an  infinitesimal  volume  of  material  in  this 
situation  is  given  by  equation  (1)  where 


AE^  = AE^  - AE^ 
DEL 


(2) 


^F.  //.  Attix,  W,  C-  Roeschf  and  E.  Tochilin,  Radiation  Dosimetry, 
2nd  Ed.,  Academic  Press,  New  York  (1966). 

^F.  H.  Attix,  Basic  Gamma-Ray  Dosimetry,  Health  Physics,  ^ (1968), 
49. 

^Robley  D.  Evans,  The  Atomic  Nucleus,  McGraw-Hill  Book  Co.,  New  York 
(1955),  723. 
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in  which  AE  is  the  energy  entering  the  mass  Am;  AE  is  the  energy 
E L 

leaving  it.  This  situation  is  pictured  in  figure  1 (a) . Imagine  a given 
(not  necessarily  microscopic)  volume  of  material  exposed  to  high-energy 
photon  radiation  and  the  need  to  know  the  dose  at  a point  P.  If  Am  is 
the  infinitesimal  mass  of  material  surrounding  P,  then  the  dose  is  given 
by  equations  (1)  and  (2)  where  equation  (2)  becomes 


AE^  = AE^(e)  - AE  (e)  + AE^(y)  - AE  (y) 

U £j  Li  Li  L 


(3) 


(a) 


(b) 


7 

7 

y 


(c) 


Figure  1.  Dose:  (a)  definition  and  applications  to  (b)  charged 

particle  equilibrium  and  (c)  Bragg-Gray  cavity  theory. 
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The  parenthetical  modifiers  refer  to  electrons  (e)  and  photons  (Y) • This 
process  is  illustrated  in  figure  1(b).  All  these  Ae*s  are  difficult  to 
identify  and  calculate  in  practice.  However,  equation  (3)  can  be 
significantly  simplified  for  a condition  that  is  of  considerable 
practical  importance  to  dosimetry,  namely,  that  of  charged  particle 
equilibrium  (CPE) . 

2.1  Charged  Particle  Equilibrium 


little  material  in  front  of  Am  from  which  electrons  can  be  scattered 
into  Am.  Considering  points  farther  to  the  right,  a point  in  the 
material  is  reached  where  as  many  electrons  are  scattered  into  Am  as  are 
scattered  out.  This  is  the  condition  for  charged  particle  equilibrium: 


Considering  figure  1(b),  when  point  P is  very  near  the  left 

side  of  the  indicated  volume,  AE  (e)  will  be  very  small,  because  there  is 

E 


AE  (e)  = AE^  (e)  . 


E L 


(4) 


When  this  condition  is  satisfied,  equation  (3)  reduces  to 


AE^  = AE^(y)  - AE^  (y) 

U ri  Li 


(5) 


for  CPE.  But,  generally. 


(6) 


where  At  is  the  thickness  of  material  traversed,  and  y /p  is  the  mass 

en 

energy  absorption  coefficient.  For  small  y /p,  equation  (6)  becomes 

en 


(7) 


If  this  result  is  put  into  equation  (5)  , 


(8) 
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for  CPE,  and  from  equation  (1) 


AE 


D 


(9) 


for  CPE , where 


ae  (y) 

£j 


Y Aa 


the  incident  photon  energy  fluence,  and  A A is  an  element  of  area. 


As  an  important  consequence  of  equation  (9) , the  ratio  of  the 


equilibri\im  dose  in  material  1 to  the  equilibrium  dose  in  material  2 , 
for  the  same  incident  monoenergetic  photon  fluence,  is  given  by 


2 . 2 Bragg-Gray  Cavity  Theory 

Assume  now  that  the  infinitesimal  volume  in  figure  1(b)  is  not 
of  the  same  substance  as  the  surrounding  material,  but  is  a small  cavity 
filled  with  some  other  material.  Let  the  cavity  material  be  1;  call  the 
surrounding  material  2.  This  is  illustrated  in  figure  1(c).  Further, 
ass\ime  that  electronic  equilibrium  (CPE)  exists,  and  that  the  dimensions 
of  the  cavity  are  small  with  respect  to  the  range  of  the  secondary 
electrons. 


*This  is  only  true  for  monoenergetic  photons.  When  there  is  a photon 
energy  spectrum , <py  CEJ  , the  correct  expression  is 


Dl 

dJ 


(10) 


for  CPE. 


where 


dE 


is  the  differential  photon  energy  fluence  spectrum. 
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Subject  to  these  restrictions,  the  ratio  of  the  dose  absorbed 
by  the  cavity  material  to  that  absorbed  by  the  surrounding  material  is 


Di  Si 

D2  S2 


(11) 


where  is  the  effective  mass  stopping  power, 

— (MeV-cm^/gram)  , 
p dx 


for  the  electron  spectrum  generated  by  the  photon  field  in  material  1 , 
This  equation  will  be  applied  to  finding  the  dose  absorbed  by  a small 
thermoluminescent  dosimeter  (TLD)  in  a thick-walled  capsule. 


3.  TYPES  OF  DOSIMETERS 

The  problem  of  what  dosimeter  to  choose  for  use  at  AURORA  was 
circumscribed  by  the  probable  requirements  of  experimenters  doing  large 
systems  tests.  These  requirements  were 

(a)  Small  size,  to  fit  into  electronic  modules  with  a high 
component  packing  density 

(b)  Capability  of  making  a large  number  of  measurements  per  shot, 
possibly  200 

(c)  Mea' 'irements  of  dose  in  rads  (Si)  or  rads(Ta) 

(d)  Useful  dose  range  from  1 to  10^  rads  (Si). 

Clearly,  the  dosimeters  had  to  be  passive,  since  instrximenting  2 00 
active  channels  with  cables  and  recording  equipment  would  have  been 
prohibitive.  The  size  requirement  ruled  out  chemical  dosimeters  and 
ionization  chambers  (not  seriously  considered,  anyway).  Finally,  the 
extremely  large  dose  range  eliminated  glass  rods,  radiachromic  dye 
films,  and  photographic  film,  leaving  the  TLD*s  as  the  only  reasonable 
choice.  They  are  a good  choice  for  other  reasons  as  well.  They  are 
available  from  several  manufacturers  in  a large  variety  of  compositions 
and  physical  form.  The  TLD  readout  instruments  are  also  available 
from  a number  of  manufacturers  as  standard  catalog  items.  The  TLD*s  are 
inexpensive  and  easy  to  use  and  retain  dose  information  for  years  at 
room  temperature.  Further,  the  quality  control  in  their  manufacture  is 
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now  such  that  dose  measurements  can  be  made  with  a precision  of  ±2  to  4 
percent  (at  1 o)  without  use  of  any  unusual  techniques  (such  as  mass 
screening  or  individual  precalibration) . 

Thermoluminescent  (TL)  phosphors  are  normally  salts  of  alkali  or 
alkaline-earth  metals,  quite  often  doped  with  trace  quantities  of  other 
elements.  Some  of  the  more  commonly  used  phosphors  are  lithium 
fluoride,  LiF;  manganese-doped  lithium  tetraborate,  Li2B407:Mn; 
manganese -doped  calcium  fluoride,  CaF2:Mn;  dysprosium-doped  calcium 
fluoride,  CaF2:Dy;  and  manganese-doped  calcium  sulfate,  CaSO^iMn.  The 
mechanism  by  which  these  materials  store  energy  proportional  to  dose  and 
release  it  later  for  readout  is  similar  for  all  of  them, 

Thermoluminescent  phosphors  are  insulators  with  electron  trapping 
centers  in  the  band  gap.  Ionizing  radiation  imparts  enough  energy  to 
valence  electrons  to  raise  them  to  the  conduction  band,  whence  many  are 
trapped.  In  a good  TL  phosphor,  these  trapped  electrons  are  stable  at 
room  temperature  and  remain  trapped  for  a long  time.  The  density  of 
such  trapped  electrons  is  proportional  to  the  energy  deposited  in  the 
material  per  unit  mass,  that  is,  the  dose.  To  read  the  dosimeter,  the 
material  is  heated  to  a sufficiently  high  temperature  to  thermally 
excite  the  trapped  electrons  into  the  conduction  band,  whence  they 
return  to  the  valence  band  with  the  emission  of  light.  This  process  is 
diagrammed  in  figure  2. 

The  light  output  as  a function  of  phosphor  temperature  is  shown  for 
a representative  phosphor  in  figure  2(b).  Several  peaks  are  seen 
corresponding  to  various  trap  depths.  Either  the  peak  intensity  or  the 
integral  under  the  curve  can  be  used  as  a measure  of  the  radiation  dose. 
When  the  integral  is  used,  low  temperature  peaks  corresponding  to 
shallow  traps  can  introduce  variations  into  the  data  by  becoming  partly 
de-excited  at  room  temperatures.  When  this  is  the  case  (as  it  is  for 
LiF) , better  reproducibility  can  be  achieved  either  by  removing  the 
peaks  by  a low  temperature  preannealing  step  or  by  electronically 
integrating  only  between  two  preset  temperatures,  Ti  and  T2 , as 
indicated  by  the  shading  in  figure  2 (b) . 

Of  all  the  phosphors  mentioned  above,  CaF2:Mn  was  chosen  for  routine 
dosimetry  at  AURORA  for  the  following  reasons.  Figure  3 shows  the 
thermoluminescence  (TL)  as  a function  of  exposure  for  the  three 
phosphors  LiF,  CaF2:Mn,  and  Li2B407:Mn.  The  CaF2:Mn  and  Li2B40y:Mn  are 
usable  at  much  higher  exposures  than  LiF , about  10^  versus  about 
2 X 10^  R for  LiF.  Also,  the  TL  response  of  LiF  departs  farthest  from 
linearity,  whereas  CaF2:Mn  is  the  most  linear.  It  has  also  been 
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Figure  2.  The  radio-thermoluminescence  process:  (a)  band 

diagram  and  (b)  typical  glow  curve. 

reported  that  the  TL  response  versus  exposure  of  LiF  depends  on  the 
exposure  rate,^  whereas  the  response  of  CaF2:Mn  and  Li2B407:Mn  is 
independent  of  rate.^ 


H.  Attix,  Dosimetry  for  the  Aurora  Facility , Final  Report ^ Naval 
Research  Laboratory,  Washington,  DC,  NRL  Memorandum  Report  2552,  Revised 
(March  1973) . 

* N.  Goldstein , Dose-Rate  Dependence  of  Lithium  Fluoride  for  Exposures 
above  15,000  R per  Pulse,  Health  Physics,  22  (1972),  90, 
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Figure  3.  Thermoluminescence  response  versus  exposure  of 
three  phosphors. 


Most  experimenters  at  AURORA  are  interested  in  having  the  dose 
measured  in  units  of  rads (Si).  The  dose  in  silicon,  with  the  dose  known 
in  some  material,  X,  is  given  by 


D 


Si 


(12) 


where  u (E) /p  is  the  mass  energy  absorption  coefficient  for  the 
en 

indicated  material  at  photon  energy  E.  If  the  photons  are  not 
monoenergetic,  and  if  the  photon  energy  spectrum  d({)(E)/dE,  [MeV/MeV-sr], 
is  known,  the  conversion  becomes 
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If  the  spectrxun,  4<(E)  , is  known  only  approximately,  as  is  usually  true, 
the  dose  conversion  can  be  calculated  only  approximately. 

The  mass  energy  absorption  coefficients  for  Si  and  CaF2:Mn  are 
almost  identical  over  the  range  300  keV  < E < 10  MeV  (fig.  4).  The 
ratio  of  equation  (12)  is  constant  at  0.99  over  the  same  range. 
Equation  (12)  then  becomes 


D^.  s 0.99  ^ .. 

Si  CaF2:Mn 


(14) 


and  this  relationship  is  independent  of  photon  energy  over  the  given 
energy  range.  The  practical  consequence  of  this  is  that  CaF2 :Mn 
dosimeters  can  be  used  to  measure  silicon  dose  with  no  need  for  detailed 
knowledge  of  the  photon  energy  or  spectrum. 

Finally,  a secondary  advantage  of  CaF2:Mn  is  that  its 
thermoluminescence  response  has  only  one  relatively  high  temperature 
glow  peak,  so  that  no  preread  annealing  steps  are  required  to  achieve 
reproducible  results. 
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Nonlinearity  of  CaFp ;Mn> — Earlier  references  to  CaF2 :Mn  have 
generally  reported  its  thermoluminescence  response  to  be  linear  with 
exposure  froin  1 R to  about  2 x 10^  R,  ^ However,  in  the  process  of 
choosing  a dosimeter  for  use  at  AURORA,  it  was  discovered  that  CaF2 :Mn 
is  indeed  not  linear  for  exposures  above  10^  R (fig.  5) . The 
differences  in  the  nonlinearity  curves  for  the  various  forms  of  CaF2 :Mn 
dosimeters  are  attributed  to  a combination  of  an  actual  supralinearity 
of  light  output  that  is  the  same  for  all  three  dosimeters,  with  a dose 
dependent  darkening  of  the  material  causing  varying  degrees  of  light 
sel f -absorption . 

If  the  TL  response  of  CaF2:Mn  had  been  truly  a linear  function  of 
dose,  it  would  have  been  possible  to  adjust  a TL  analyzer  so  that  its 
output  (i.e.,  light-sum  response)  was  numerically  equal  to  the  dose. 
The  consequences  of  having  to  ^pply  a dose— dependent  nonlinearity 
correction  to  all  readings  are  discussed  fully  in  section  6 on  data 
reduction . 


Figure  5.  Thermoluminescence  light-sum  response  per  roentgen 
for  CaF2 :Mn  (from  S.  G.  Gorbics,  F.  H.  Attix,  and 
K.  Kerris,  Health  Physics,  25  (1973),  499). 


Fading  of  CaF2 :Mn. — A second  property  of  CaF2:Mn  that  has  an 
important  effect  on  data  analysis  is  that  the  light  output  for  a given 
dose  is  not  independent  of  the  elapsed  time  between  irradiation  and 
readout.  This  fading  of  the  TL  response  is  shown  in  figure  6. 


3 

F.  Hm  Attix f W.  C,  Roesch , and  E,  Tochilin , Radiation  Dosimetry , 
2nd  Ed.,  Academic  Press,  New  York  (1966). 
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TIME  AFTER  IRRADIATION  (hr) 

Figure  6.  Fading  of  the  thermoluminescence  response  of  CaF2:Mn. 


If  Yq  is  the  light-sum  response  of  the  TLD  when  read  at  time  tg 
after  irradiation,  then  the  light- sum  response  at  other  times  is  quite 
well  represented  by^/* 


y = Yo  (tAo)' 


.-0.017 


(15) 


Since  the  elapsed  time  t varies  over  a number  of  hours  or  days  when  a 
large  number  of  TLD's  is  exposed  on  a single  shot,  this  fading 
correction  cannot  be  ignored,  but  must  be  included  in  the  data. 

CaF2:Mn  TL  phosphor  is  commercially  availcible  in  three  forms: 

(a)  80-200  mesh  powder^  I 

(b)  0.125  X 0.125  X 0.035  in.  hot-pressed  blocks  I 

(c)  6-mm  diam,  0. 4-mm-thick  discs  of  Teflon  intimately  mixed 
with  finely  powdered  CaF2:Mn.’t* 


G.  GorbicSr  F.  H.  Attix,  and  K.  Kerris,  Thermoluminescent 
Dosimeters  for  High-Dose  Applications,  Health  Physics,  2^  (1973),  499. 

'^Reference  2 gives  y/y^  = 1 - 0.017Zn^/t^  Equation  (15)  represents 
the  data  equally  well. 

’^Teledyne  Isotopes,  Inc.,  Westwood,  NJ . 

^-Harshaw  Chemical  Co.,  Solon,  OH. 
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The  decision  on  which  one  to  use  for  routine  dosimetry  at  AURORA  was 
based  principally  on  matters  of  convenience  and  cost.  The  detailed 
tradeoffs  ^ are  summarized  below. 

Powder  had  the  lowest  price  per  dosimeter  (10<?)  , assuming  10-mg 
samples  and  no  reuse  of  the  powder,  but  its  use  was  considered 
impractical  because  of  the  time  involved  in  measuring  out  and 
encapsulating  uniform  10-mg  samples  of  powder.  It  is  also  very 
difficult  to  envision  an  automatic  sample  changer  for  a TLD  reader  that 
could  hcindle  powder  samples. 

Hot-pressed  blocks  had  the  highest  unit  price  ($1.75), 
requiring  their  being  annealed  and  reused  repeatedly.  They  can  be  if 
they  are  never  irradiated  above  50  to  100  rads (Si).  However,  at  higher 
doses,  their  calibration  changes  after  several  irradiation  and  annealing 
cycles,  making  this  mode  of  operation  impractical . 

The  Teflon  discs  were  priced  low  enough  (30<^)  to  permit  their 
being  used  once  only.  They  are  easy  to  handle  and  adaptable  to 
automated  loading.  Finally,  the  sample-to-sample  reproducibility  for 
TLD's  from  the  same  production  lot  was  ±3.1  percent  (1  a)  , which  is 
acceptable  for  routine  use  at  AURORA.  The  scimple- to-sample 
reproducibility  gets  quite  a bit  worse  around  10  rads  (Si),  so  that  these 
TLD's  are  not  used  for  measuring  doses  below  50  rads  (Si). 

The  following  TLD's  were  adopted  as  standard  for  routine  use  at 

AURORA: 


(a)  For  50  D 10^  rads  (Si)  : Teflon  disc  containing  5 
percent  by  weight  CaF2:Mn;  Teledyne  Isotopes  SD-CaF2 :Mn-0.4L.  These 
TLD's  are  procured  in  quantities  of  15,000  from  the  same  production  lot. 

(b)  For  mrad(Si)  D 50  rads(Si)  : Harshaw  Chemical  Company, 
0.125  X 0.125  X 0.035  in.  hot-pressed  CaF2:Mn  blocks.  One  thousand  were 
procured  from  the  same  production  lot.  It  is  not  known  how  often  these 
can  be  irradiated,  annealed,  and  reused  before  the  individual 
calibrations  drift  enough  to  make  a batch  calibration  invalid.  At  very 
low  doses  (D  1 rad (Si)) , dry  nitrogen  gas  flow  must  be  used  with  these 
TLD*s  during  the  readout  cycle  to  suppress  spurious  luminescence  caused 
by  absorbed  oxygen. 


^ F.  H.  AttiXf  Dosimetry  for  the  Aurora  Facility , Final  Report^  Naval 
Research  Laboratory,  Washington,  DC,  NRL  Memorandum  Report  2552,  Revised 
(March  1973) . 


17 


4.  TYPES  OF  MEASUREMENT 


These  TLD*s  can  be  used  to  make  two  entirely  different  kinds  of 
silicon  dose  measurements,  and  it  is  very  important  to  understand  the 
distinction  between  them. 

Assume  that  an  experimenter  is  irradiating  a large,  complex 
electronic  system  at  AURORA.  Somewhere  inside  this  system  in  some 
critical  circuit  is  a silicon  transistor  or  integrated  circuit,  and  the 
experimenter  wants  to  know  the  dose  absorbed  by  this  component  in 
rads (Si).  The  active  volume  of  silicon  is  very  small  compared  to  the 
rest  of  the  system,  so  that  all  the  energy  deposited  in  this  volume  of 
silicon  comes  from  secondary  electrons  produced  by  photons  somewhere  in 
the  surrounding  hardware,  and  none  comes  from  primary  photon 
interactions  in  the  silicon  itself. 

The  way  to  determine  this  dose  is  to  put  a small  dosimeter  very  near 
this  critical  component.  It  will  then  see  very  nearly  the  same 
secondary  electron  environment  as  the  transistor,  and  the  dose  absorbed 
by  the  TLD  will  be  related  to  the  transistor  dose  by  the  ratio  of  the 
two  electron  stopping  powers.  Since  the  dosimeter  has  a stopping  power 
almost  identical  to  that  of  silicon  over  a wide  range  of  electron 
energies,  this  ratio  is  nearly  one,  and  it  is  constant  over  a large 
energy  range,  making  the  measurement  largely  independent  of  any  detailed 
knowledge  of  the  AURORA  bremsstrahlung  spectnam. 

So  to  measure  the  local  silicon  dose  inside  a nonuniform,  rather 
lumpy  system,  a bare  CaFarMn  TLD  is  placed  as  closely  as  possible  to  the 
exact  point  of  interest.  This  kind  of  measurement  is  normally  called 
Type -A  at  AURORA. 

However,  assiame  that  this  is  not  the  measurement  the  experimenter 
wants.  Rather,  he  wants  to  know  the  photon  fluence  incident  on  his 
system  in  rads  (Si)  . Of  course,  photon  fluence  is  not  measxired  in 
rads  (Si)  but  MeV/cm^  as  a function  of  photon  energy.  So  what  does 
rads (Si)  incident  mean?  It  means,  simply,  the  dose  that  would  be 
deposited  in  silicon  under  conditions  of  CPE  by  the  incident  fluence, 
whatever  it  may  be.  The  actual  fluence  as  a function  of  photon  energy 
would  be  a very  difficult  measurement,  but  the  equilibrium  silicon  dose 
is  easy.  One  simply  surrounds  a silicon  dosimeter  with  enough  silicon 
to  produce  CPE.  In  actual  practice,  one  uses  a CaF2:Mn  TLD  and 
surrounds  it  with  aliaminiam.  Doing  so  requires  making  a photon 
absorption  correction  for  the  aluminum  and  a stopping  power  correction 
for  the  CaF2,  but  these  corrections  are  very  nearly  unity  and  constant 
for  300  keV  ^ hv  ^10  MeV,  so  that  this  procedure  is  completely  valid 
and  independent  of  any  detailed  knowledge  of  the  AURORA  spectrum. 
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How  much  aluminum  is  required  to  guarantee  CPE?  This  question  was 
answered  experimentally  by  determining  the  dose  absorbed  by  the  TLD  as  a 
function  of  aluminum  thickness^  (fig.  7).  The  initial,  rather  sharp 
rise  shows  that  in  this  region,  more  electrons  are  scattered  out  of  the 
center  of  the  aluminum  capsule  than  are  scattered  into  it  from  the 
walls.  The  secondary  electron  flux  is  low  because  there  is  no  CPE. 
Charged  particle  equilibrium  exists  for  all  points  to  the  right  of  the 
peak,  but  the  dose  slowly  decreases  because  of  attenuation  of  the 
photons  by  the  aluminum.  The  optimum  aluminum  thickness  to  assure  CPE 
without  introducing  attenuation  is  seen  to  be  1 gram/cm or  0.146  in. 

■V 

So,  to  measure  the  silicon  equilibrium  dose  in  the  free  field,  one 
surrounds  the  TLD  with  an  aluminum  capsule  of  5/32  in.  (0.156  in.)  wall 
thickness.  This  type  of  measurement  is  normally  referred  to  at  AURORA 
as  Type-B.  Figure  8 is  an  engineering  sketch  of  the  standard  AURORA 
Type-B  capsule. 

A third  kind  of  measurement  that  is  sometimes  requested  by  AURORA 
experimenters  is  the  tantalum  equilibrium  dose,  rads(Ta).  This  dose  is 
measured  by  a small  tantalum  dosimeter  surrounded  by  an  equilibrium 

ALUMINUM  THICKNESS  (grom/cm^) 


0 1 0 2 0.3  0.4 

ALUMINUM  THICKNESS  (in.) 


Figure  7.  Aluminum  depth  dose  (F.  H.  Attix,  NRL  Memorandum 
Report  2552  (March  1973)). 


^F.  H.  Attix,  Dosimetry  for  the  Aurora  Facility,  Final  Report,  Naval 
Research  Laboratory , Washington , DC,  NRL  Memorandum  Report  2552,  Revised 
(March  1973) . 
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Figure  8.  Standard  Type-B  dosimeter  capsule - 


thickness  of  tantalum.  Unfortunately/  at  the  present  time,  there  are  no 
high-Z  TLD  phosphors  that  approximate  the  electron  stopping  power  of 
tantalum.  In  principle,  it  would  be  possible  to  use  a CaP2  TLD  with  an 
equilibrium  shield  of  tantalum  and  to  make  a stopping  power  correction 
for  the  CaF2.  Unfortunately,  the  CaF2  TLD  is  not  sufficiently  small 
compared  with  the  tantalum  shield  for  simple  Bragg -Gray  cavity 
theory^  to  apply.  What  is  done  at  AURORA  is  to  calibrate  the 
tantalum-clad  TLD's  against  a tantalum  calorimeter. 

The  tantalum  capsules  are  formed  from  a 0.015-in. -thick  tantalum 
sheet.  This  configuration  is  called  Type-C.  Figure  9 (p.  21)  shows  both 
the  Teflon  disc  CaF2:Mn  TLD's  and  the  solid  CaF2:Mn  blocks  in  their 
Types  A,  B,  and  C configurations.  Type-C  TLD  capsules  are  7/16  in. 
square  by  1/16  in.  thick. 

5 . INSTRUMENTATION 

5.1  Dose  Measurement 

Apparatus  for  reading  out  TLD's  typically  consist  of  a heating 
element,  circuitry  that  provides  a reproducible  temperature  cycle,  and  a 
photomultiplier  tube  and  amplifiers  to  detect  the  light  emitted  during 

^F.  H*  Attix,  W.  C.  Roescb^  and  F.  Tochilinr  Radiation  Dosimetry , 
2nd  Ed.,  Academic  Press,  New  York  (1966). 
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DOSIMETERS  SOLID  Co Mn  DOSIMETERS 


TYPE -A  CONFIGURATION 

(negative  No.  099-74) 

Figure  9.  Types  A,  B,  and  C dosimeter  configurations. 

the  temperature  cycle.  The  heating  element  is  usually  a small,  shallow 
dish  made  of  some  suitable  material  such  as  Monel  Metal,  Nichrome , or 
platinum,  which  is  heated  when  current  is  passed  through  it. 
Temperature  cycling  can  be  accomplished  by  a programmed  power  supply, 
often  with  temperature  feedback  from  a thermocouple  in  contact  with  the 
heating  element. 

The  light  detecting  circuitry  generally  includes  a filter  to 
block  infrared  radiation,  a visible  light  filter  to  select  the  proper 
spectral  band  for  the  phosphor  in  use,  a photomultiplier  tube,  and  a 
stcible  high  voltage  power  supply.  One  most  generally  determines  the 
sensitivity  of  the  reader  by  adjusting  the  photomultiplier  voltage.  The 
photomultiplier  current  can  be  amplified  and  recorded  on  a strip-chart 
recorder  to  show  the  glow  curve  or  integrated  electronically.  Either 
the  peak  of  the  glow  curve  or  the  integral  under  the  curve  can  be  used 
as  a measure  of  the  dose  absorbed  by  the  TL  material. 

The  primary  criteria  for  the  design  of  the  readout 
instrumentation  to  be  used  at  AURORA  were  operational  simplicity  and  a 
high  degree  of  automation.  In  particular,  all  steps  requiring  the 
operator  to  copy  data  manually  were  to  be  eliminated,  with  everything 
being  designed  for  eventual  fully  automatic  hands-off  operation.  These 
criteria  were  imposed  because  of  the  large  volume  of  dosimetry  readings 
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required  by  many  AURORA  users.  During  a large  systems  test,  which  might 
last  for  several  weeks,  one  person  is  kept  busy  all  day  just  reading 
TLD's.  It  was  therefore  considered  quite  important  to  design  the 
equipment  to  reduce  operator  fatigue  as  much  as  possible  and  to 
eliminate  those  steps  in  which  operator  fatigue  might  introduce  errors 
into  the  data. 

A secondary  requirement  on  the  TLD  readout  instrumentation  was 
that  it  be  flexible  enough  to  accommodate  the  large  variety  of  TL 
materials  required  for  specialized  measurements.  Although  CaF2:Mn  is 
the  most  commonly  used  phosphor  at  AURORA,  it  was  thought  possible  that 
LiF,  CaS04,  Li2B407,  and  other  phosphors  would  be  used  for  some 
applications. 

The  instrument  chosen  for  use  at  AURORA  was  the  Harshaw  Model 
2000  Thermoluminescence  Analyzer.  It  has  a very  versatile  temperature 
cycle  programmer  cind  linear  and  logarithmic  glow- curve  outputs;  these 
fully  satisfy  the  flexibility  criterion.  The  integral  under  the  glow 
curve  (light  sum)  is  displayed  digitally  and  is  also  available  in 
digital  form  at  a rear  panel  connector,  making  it  possible  to  achieve 
semiautomated  or  fully  automatic  operation. 

The  most  likely  chcince  for  operator  error  in  making  many 
repetitive  measurements  is  in  miscopying  the  data.  Therefore  the  first 
step  toward  automating  the  readout  process  was  to  add  a printer  to  print 
the  light  sum  of  each  measurement  automatically.  The  Harshaw  TLD  reader 
furnishes  a print  commcind  at  its  rear  panel  connector  at  the  end  of  each 
readout  cycle,  so  that  the  printer  connection  is  easily  made.  The 
printer  used  for  this  purpose  is  a Systron-Donner  Model  5103  20-column 
printer. 


The  sensitivity  of  the  TLD  reader  is  adjusted  so  that  the  value 
of  the  light  sum  as  printed  on  the  paper  tape  is  approximately  equal  to 
the  dose  in  rads (Si),  except  for  a factor  equal  to  a power  of  10.  Since 
the  dose  response  of  CaF2:Mn  is  not  exactly  linear,  the  printer  output 
can  be  made  equal  to  the  dose  only  at  one  particular  dose.  In  practice, 
the  calibration  is  chosen  to  give  a one-to-one  correspondence  at  a dose 
of  about  20  krads(Si),  because  TLD*s  are  most  frequently  used  to  measure 
doses  from  5 to  40  krads(Si)  at  AURORA. 

Another  difficulty  was  finding  a way  to  print  the  proper 
multiplier  (i.e.,  power  of  10)  for  the  printed  data.  On  its  lowest 
range,  the  TLD  reader  rear  panel  output  consists  of  four  digits,  with 
the  least  significant  digit  equal  to  1 pC.  The  display  automatically 
ranges  up  by  factors  of  100;  that  is,  when  the  display  reaches  9999  pC, 
the  addition  of  one  more  picocoulomb  results  in  0100  (xlOO)  pC,  and  so 
on  up  through  three  rcinge  changes.  The  proper  range  is  automatically 
indicated  on  the  front  panel  by  a shifting  decimal  point  and  illuminated 
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nanocouloitib  and  microcoulomb  indicators.  However,  at  the  rear  panel 
connector,  the  range  in  use  is  indicated  by  a two-digit  binary  number, 
OO2  signifying  the  lowest  range  (1  digit  = 1 pC)  and  II2  indicating  the 
highest  (1  digit  = 1 pC) . Further,  the  basic  calibration  of  the  reader 
can  be  1,  10,  or  100  rads (Si) /nC. 

This  problem  was  solved  by  constructing  an  interface  circuit 
that  does  two  things.  First,  it  puts  a constant  decimal  point  between 
the  most  significant  and  next  most  significant  digits  printed  on  the 
paper  tape,  thus:  X.XXX.  Secondly,  it  accepts  the  two-digit  range  code 
from  the  Harshaw  reader  and  a two-digit  code  (programmed  by  toggle 
switches)  indicating  the  basic  reader  sensitivity.  It  decodes  this 
input  and  commands  the  printer  to  print  a fifth  digit  representing  the 
correct  exponent  for  the  power-of-10  multiple  to  be  applied  to  the  data. 
For  example,  a printout  of  1.234  5 R is  to  be  interpreted  as 
1.234  X 10^  rads(Si),  except  for  nonlinearity  and  fading  corrections. 

The  truth  table  for  these  decoding  operations  is  given  in 
table  I.  An  example  may  illustrate  how  this  wor)cs.  Suppose  a reader 
with  a basic  calibration  of  10  rads (Si) /nC  indicates  a light  sum  of 
03.22  pC  on  its  front  panel  display.  The  printer  output  will  be  0322 
and  the  range  indication,  IO2.  For  a toggle  switch  setting  of  OI2,  the 
decoded  exponent  of  10  will  be  5,  and  the  printer  will  print  0.322  5 R 
(0.322  X 10  5 rads  (Si))  .' 

(Verify:  10  rads(Si)/nC  x 3.22  x 10^  nC  = 3.22  x 10**  rads(Si)). 


TABLE  I.  RANGE  DECODING  TRUTH  TABLE 


a 

Harshaw  range  data 

Reader 

OO2 

(1  rad(Si)/nC)  (10 

calibration  toggle 
OI2 

rads (Si) /nC)  (100 

switch  setting 
IO2 

rads(Si)/nC)  (1000 

112 

rads (Si)/nC) 

OO2 

(LSD^  = 1 pC) 

0 

Decoded  exponent 
1 2 

3 

OI2 

(LSD  = 100  pC) 

2 

3 

4 

5 

10  2 

(LSD  = 10  nC) 

4 

5 

6 

7 

II2 

(LSD  = 1 UC) 

6 

7 

(c) 

- 

^Subscript  2 * betse  - 2 (binar\^)  numbers. 
^LSD  = least  significant  digit. 

^ = no  code. 
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The  decoding  circuit  is  easily  implemented  with  an  SN74154 
4“  to  16-line  decoder  and  three  SN7430  eight-input  NAND  gates  (fig.  10) . 
Two  toggle  switches  furnish  either  a logic  0 or  1 at  the  A and  B inputs. 
The  Harshaw  range  output  is  connected  to  the  C and  D inputs.  The  output 
to  the  printer  is  binary-coded  decimal  (BCD) , with  the  most  significant 
bit  hardwired  to  logic  0.  The  BCD  output  of  the  circuit  of  figure  10 
follows  the  truth  table  of  table  I. 

To  identify  the  numbers  on  a long  strip  of  printer  tape  with 
specific  dosimeters,  a provision  to  print  a three-digit  identifying 
number  along  with  each  dose  datum  was.  incorporated  into  the  interface 
unit.  This  identifying  number  is  input  via  a 10-key  keyboard  and 
printed  at  the  extreme  left  of  the  paper  strip;  so  a datum  of 


357  0.322  5 R 


means  that  the  dose  (uncorrected)  of  3.22  x lo^  rads (Si)  is  to  be 
identified  with  TLD  number  357.  The  identifying  number  can  be  keyed  in 
two  modes,  selected  by  a toggle  switch.  In  the  MANUAL  mode,  used  when 
there  are  few  TLD*s  or  the  identification  numbers  are  in  some  random 
order,  the  number  register  is  cleared  after  each  print  cycle,  and  a new 
number  must  be  keyed  in  for  each  reading.  If  there  are  many  readings  in 


SN  7430 

8- INPUT  NANO  GATES 


Figure  10.  Range  decoding  circuit. 
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numerical  sequence,  an  AUTOMATIC  mode  is  provided.  In  this  mode,  only 
the  first  number  need  be  keyed  in;  the  identifying  number  is 
automatically  increased  by  one  every  time  a reading  is  made  and  another 
READ  cycle  is  initiated.  This  increase  is  made  possible  by  a READ  level 
change  being  available  at  the  rear  panel  connector  of  the  Harshaw 
reader.  A light-emitting  diode  display  on  the  interface  unit  displays 
the  identifying  number  and  the  exponent  of  10. 

The  TLD  identifying  number  circuit  is  shown  in  figure  11. 
Whenever  a digit  key  is  pressed,  the  BCD  data  appear  at  the  data  input 
of  all  three  SN74192  asynchronous  BCD  up-down  counters.  If  this  is  the 
first  key  that  is  pressed,  the  address  shift  register  applies  a load 
enable  pulse  to  the  counter  for  the  most  significant  digit  (MSD) . The 
second  time  a key  is  pressed,  the  second  digit  is  enabled,  and  the  third 
keystroke  enables  the  least  significant  digit  (LSD) . The  outputs  of  the 
counters  are  connected  to  the  first  three  columns  of  the  printer.  Upon 
completion  of  a PRINT  cycle,  the  digital  printer  sends  a pulse  to  the 
mode  control  circuitry.  If  the  switch  is  in  the  MANUAL  position,  this 
pulse  clears  the  BCD  counters  and  resets  the  address  shift  register.  If 
the  switch  is  in  the  AUTOMATIC  position,  the  RESET  pulse  is  blocked,  and 
the  next  succeeding  READ  cycle  initiation  causes  a COUNT  pulse  to  be 
applied  to  the  LSD  counter. 


MSD  LSD 


AUTOMATIC 

MODE 

CONTROL 


Figure  11.  Dosimeter  identifying  number  circuit. 
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The  address  shift  register  consists  of  two  SN7474  dual-D 
flip-flops,  one  SN7400  quad  two-input  NAND  gate,  and  two  SN74121 
one-shots.  The  mode  control  circuit  contains  two  SN74121  one-shots. 
The  entire  interface  unit  circuit  is  shown  in  figure  12.  LOAD  ENABLE 
pulses  are  triggered  by  the  leading  edge  of  the  keyboard  STROBE  pulse; 
the  SHIFT  REGISTER  CLOCK  pulse  is  derived  from  the  trailing  edge  of  this 
STROBE  pulse.  RESET  and  CLEAR  pulses  are  generated  by  the  trailing  edge 
of  either  the  keyboard  CLEAR  KEY  pulse  or  the  printer  +INHIBIT  pulse. 
The  COUNT  UP  one-shot  puts  out  a 1-s  pulse  to  suppress  multiple  pulsing 
by  a very  bouncy  READ  switch  on  the  Harshaw  reader.  The  mode  is 
controlled  by  switching  the  logic  levels  of  appropriate  inputs  on  the 
RESET  and  COUNT  UP  one-shots. 


Figure  13  shows  the  entire  TLD  reader  system.  The  only  manual 
operations  required  are  opening  the  sample  drawer,  inserting  a TLD, 
closing  the  sample  drawer,  and  pressing  the  READ  button.  Data  recording 
is  automatic,  eliminating  the  possibility  of  erroneous  data  entry  by 
hand.  A portion  of  a sample  data  tape  is  shown  below. 


0 1 8 

0 4 

14  5 9 

1 

.6  6 6 

4 

R 

0 1 7 

0 4 

14  5 7 

1 

.5  9 6 

4 

R 

0 1 6 

04 

14  4 5 

1 

.5  0 5 

4 

R 

0 1 5 

0 4 

14  4 4 

1 

.6  6 4 

4 

R 

0 1 4 

04 

14  4 3 

1 

.5  9 9 

4 

R 

0 1 3 

04 

14  4 1 

1 

.7  1 5 

4 

R 

0 1 2 

04 

14  4 0 

1 

.4  9 7 

4 

R 

0 1 1 

04 

14  3 9 

1 

.6  3 4 

4 

R 

is  the 

TLD 

identifying  number. 

The 

the  month 

and  24- 

-hr 

time 

of 

th 

The  first  number  is  the  TLD  identifying  number.  The  next  two  numbers 
are  the  day  of  the  month  and  24-hr  time  of  the  reading.  This 
information  is  necessary  for  calculation  of  the  fading  correction  for 
CaF2:Mn  TLD*s.  This  information  is  generated  by  an  Option  D102-H  clock 
board  inside  each  printer.  The  final  numbers  represent  the  dose  in 
approximate  rads (Si) . 


It  is  tempting  to  imagine  a fully  automatic  system,  including 
some  sort  of  mechanical  TLD  changing  device  that  would  eliminate  the 
steps  of  drawer  opening  and  closing  and  button  pushing.  Such  a device 
is  now  under  construction.  When  complete,  it  will  accept  a magazine  of 
up  to  200  TLD*s,  read  them,  and  record  the  data,  all  automatically. 
Since  each  complete  readout  cycle  takes  90  s,  each  automatic  readout  of 
one  200-TLD  magazine  would  represent  a savings  of  5 man-hours. 
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Figure  12,  Interface  unit  circuit 


(negative  No.  099-74) 

Figure  13.  Thermoluminescent  dosimeter  reader,  printer, 
and  interface  unit. 


5 . 2 Dose  Rate  Measurement 


The  AURORA  bremsstrahlung  pulse  is  measured  on  every  shot  by 
use  of  a scintillator-photodiode  combination.  The  detector  is  usually 
placed  so  that  the  scintillator  is  as  near  the  experimenter's  cassette 
and  the  machine  center  line  as  feasible.  The  light  from  the 
scintillator  is  coupled  to  the  photodiode  by  a hollow  aluminum  light 
pipe  4-  or  6-ft  long.  The  photodiode  is  usually  placed  near  the  floor 
of  the  test  cell,  inside  a 4-in.  thick  lead  shield.  This  arrangement  is 
diagrammed  in  figure  14. 

The  scintillator  is  usually  Pilot-B,  either  a cylinder 
3 in,  diam  x 3 in.  long  or  a disc  3 in.  diam  x 1/16  in.  thick,  depending 
on  the  dose  rate.  When  it  is  desired  to  look  at  very  high  dose  rates 
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BREMSSTRAHLUNG 


Figure  14-  Scintillator-photodiode  detector  setup- 

(>5  X 10^^  rads (Si) /s),  the  Cerenkov  radiation  from  a cylinder  of  lucite 
3 in.  dieun  x 3 in.  long  is  detected. 

The  photodiode  detector  is  powered  by  a 5-kV  power  supply,  and 
its  output  signal,  properly  attenuated,  is  displayed  by  one  or  more 
Tektronix  519  oscilloscopes.  The  circuit  for  the  photodiode  detector  is 
shown  in  figure  15.® 


VARIABLE  son  5on*  i2sn 
ATTENUATOR  MATCHING  PAD 


Figure  15.  Detector  circuit. 


^J.  F.  Doyle  and  W.  C.  Dickinson,  The  Cerenkov  Response  of  Lucite  and 
Quartz  to  Gamma  and  Fast  Neutron  Radiation,  University  of  California 
Lawrence  Radiation  Laboratory,  UCRL-7032  (TID4500)  (July  1962). 

K.  Durkee,  Design  and  Test  of  a Modular  Scintillation  Detector 
System,  EG&G , Santa  Barbara  Division,  Technical  Report  S-466-R  (March 
1970) . 


The  detector  is  not  normally  calibrated.  The  dose  rate  at  each 
spatial  point  of  interest  is  determined  by  measuring  the  dose  at  that 
point  with  a TLD,  assuming  that  the  pulse  shape  at  that  point  is  the 
Scime  as  that  at  the  Pilot--B  scintillator.  The  dose  rate  is  then 
calculated  from  the  known  dose  and  pulse  shape. 


6.  DATA  ANALYSIS 

The  number  printed  by  the  digital  printer  is  proportional  to  the 
light-sum  output  of  the  TLD.  As  such,  the  number  is  not  very  useful^ 
since  the  light  sum  (for  the  Teflon  disc  TLD*s)  is  a nonlinear  function 
of  the  dose.  In  addition,  the  light  sum  fades  with  time.  If  the 
light  sum  were  a linear  function  of  dose,  and  if  all  TLD*s  were  always 
read  at  the  same  time  relative  to  irradiation,  the  printed  number  could 
always  be  made  equal  to  the  dose  by  appropriate  adjustment  of  reader 
sensitivity  (via  the  photo-tube  high  voltage)  and  interface  unit  decoder 
switch  settings. 

As  things  stand,  however,  the  printed  number,  P,  must  be  corrected 

by  multiplying  by  a nonlinearity  correction,  K , and  a fading 

L 

correction,  K^.  The  nonlinearity  correction  is  a function  of  the  dose, 
and  the  fading  correction  is  a function  of  the  time  elapsed  between 
irradiation  and  reading. 

6.1  Nonlinearity  Correction 

We  define  the  following  quantities, 
y = light-sum  output  of  TLD 
X = exposure,  roentgens 
D . = dose,  rads (Si) 

o X 

P = TLD  reader  printout. 

We  are  looking  for  a nonlinearity  correction  factor,  K^,  such  that 

D . = K (P)  X p . (16) 

Si  L 

Figure  5 shows  the  relative  thermoluminescence,  y(X)/y(XQ),  per 
roentgen,  as  a function  of  the  Co^^  exposure,  X.  This  nonlinearity 
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curve  was  determined  by  irradiating  a number  of  TLD's  to  various 
exposures,  X,  and  determining  the  corresponding  values  of  y(X).  To  find 
K^(P) , use  the  data  in  figure  5,  but  plot  Xy(Xo)/y(X)  as  a function  of 
y(X).  Now,  P(X)  is  proportional  to  y(X),  so  P(X)  = Cy(X).  Also, 
(AURORA)  = 0.870X (Co^O) . Applying  these  two  relationships  gives  a 
new  plot  of  Dg^/P(D)  as  a function  of  P.  But  the  constant  of 

proportionality,  C,  between  the  light-sum  response  and  the  printout  is 
arbitrary,  depending  on  the  sensitivity  setting  of  the  TLD  reader. 

To  make  the  printed  data,  P,  nearly  equal  to  for  most  cases 

of  practical  interest  at  AURORA,  choose  C such  that  “ 1.000 

for  P = 1.55  X 10^*.  The  resulting  curve  of  K versus  P is  given  in 

1j 

figure  16.  This  is  our  standard  nonlinearity  correction  curve.  When 
P = 1.55  X 10^,  =1.55  X 10^  rads (Si) . All  other  readings,  P,  must 

be  multiplied  by  to  find  the  dose.  For  instance,  if 

P = 3.0  X 10 ^ for  a given  dosimeter, 

D . = K (P)  X P = 1.11  X 3 X 1q3  = 3.33  x 10^  rads(Si). 

Si  L 
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Figure  16. 
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TLD  READER  PRINTOUT,  P 

Standard  nonlinearity  correction  for  Teledyne- 
Isotopes  SD“CaF2 :Mn“0.4L. 
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This  can  be  calculated  by  hand  with  figure  16.  However,  for 

many  values  of  P,  hand  calculation  becomes  excessively  laborious.  For 

machine  calculation,  it  would  be  very  desirable  to  have  an  analytic 

expression  of  K as  a function  of  P.  The  expression  was  found  by 

leaist-squares  fitting  a fifth-degree  polynomial  to  the  data  of 

figure  16.  That  is,  represent  K by 

L 

S 

K (p)  = V'  a (log  P)"  (17) 

Li  n 

n=0 


where  the  coefficients  ao  to  as  are  found  to  be 

ao  = 0.99635 
a^  = 0.76943 

a2  = -0.67278  (18) 

as  = 0.26418 
a4  = -0.050483 
as  = 0.0036386  . 


This  polynomial  representation  deviates  from  figure  16  by  at  most 
2 percent  and  on  the  average  by  less  than  1 percent.  Given  a value,  P, 
the  silicon  dose,  is  given  by 


5 


n=0 


where  ao  to  as  are  given  by  equations  (18)  . 


(19) 


6.2  Fading  Correction 

The  analytical  expression  for  the  fading  of  the  light-sum 
response  was  given  in  equation  (15)  by 


Y = Yo(t/to)  -0-017^ 
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But  P is  proportional  to  y,  so 


P 


-0.017 


(20) 


The  reader  is  calibrated  to  give  the  correct  dose  when  the  TLD  is  read 
exactly  1 hr  after  irradiation.  So  the  P to  be  used  in  equation  (19)  is 
Pg,  which  is  related  to  the  actual  P by 


PO 


= p(t/to) 


+0.017 


(21) 


Combining  the  nonlinearity  correction  equation  (19)  with  the 
fading  correction  equation  (21)  and  setting  tQ  = 1 hr  gives  the  final 
expression  for  the  silicon  dose  absorbed  by  a TLD,  which  gives  a reading 
P when  read  at  a time  t hours  after  irradiation 


5 

D = P y'  a (log  P)"  . (22) 

ox  / ^ n 

n=0 


6.3  Solid  CaF2:Mn  TLD*s 


The  response  for  these  TLD*s  is  almost  linear  with  dose,  at 
least  over  the  range  of  0 to  200  rads  (Si).  However,  the  fading 
correction  must  still  be  made,  as  for  the  Teflon  discs.  Using  a Co^^ 
calibration  technique  exactly  as  described  for  the  Teflon  disc  TLD’s, 
and  with  the  TLD  readers  calibrated  so  that  equation  (22)  is  true,  one 
can  find  the  analogous  equation  for  solid  CaF2 :Mn  TLD's 


Si 


n=0 


b (P 
n 


^0.017jn 


(23) 


where  bo  through  b2  are  found  by  fitting  a quadratic  equation  to  the 
calibration  curve  by  the  method  of  least  squares. 

6.4  Dose  Rate  Calculation 

The  dose  rate  is  found  from  the  dose,  by  use  of  the  known  pulse 
shape  of  the  radiation  pulse  as  determined  by  the 
scintillator-photodiode  detector.  The  area  under  the  pulse  is  found  by 
numeric  integration.  Dividing  the  area  by  the  peak  height  gives  the 
effective  pulse  width.  "Effective"  pulse  width  means  the  width  of  a 
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rectangular  pulse  having  the  same  area  and  the  same  peak  height  as  has 
the  AURORA  pulse.  Since  the  dose  is  proportional  to  the  area  under  the 
AURORA  pulse,  the  dose  divided  by  the  effective  width  equals  the  dose 
rate  at  the  peak  of  the  pulse : 


dose  rate  at  peak  = 


D 


Si 


'eff 


(rads (Si)/s) 


(24) 


where 


t 


eff 


area  under  pulse 
peak  height 


(25)' 


6.5  Calculation  of  Data 


The  calculations  required  for  equations  (22)  to  (25)  appear  to 
be  quite  involved;  however,  they  are  simply  implemented  in  practice. 
All  calculations  are  done  on  the  Hewlett-Packard  9820A  programmable  desk 
calculator  equipped  with  a 9861A  typewriter,  9862A  plotter,  and  9864A 
digitizer. 

The  oscilloscope  trace  representing  the  photodiode  pulse  is 
traced  with  the  digitizer  cursor.  The  calculator  computes  the  peak 
height  and  the  area,  using  the  trapezoidal  rule,  while  the  plotter  draws 
an  enlarged  picture  of  the  pulse  on  the  data  sheet.  The  operator  then 
places  the  data  sheet  in  the  typewriter.  He  enters  the  shot  n\imber, 
time  of  irradiation,  time  of  reading  (colvimn  2 of  the  data  tape)  , and 
first  and  last  TLD  identification  numbers  (colximn  1 of  the  data  tape)  . 
He  then  enters  the  raw  data,  P (column  3 of  the  data  tape) , in  sequence. 
The  calculator  calculates  the  elapsed  time,  dose,  and  dose  rate  by 
equations  (22)  to  (25)  and  causes  the  typewriter  to  type  all  data  on  the 
data  sheet. 

Figure  17  reproduces  a typical  data  sheet.  Note  the  picture  of 
the  pulse  to  the  scale  of  100  ns/in. , pulse  parameters,  and  dose  and 
dose  rate  data.  The  raw  data  used  to  produce  this  data  sheet  are  like 
those  shown  on  the  sample  data  tape  in  section  5.  The  raw  data,  P, 
are  approximately  equal  to  the  doses,  permitting  preliminary  data 
evaluation  as  the  TLD*s  are  read. 


This  data 
programmable  desk 
data  taking.  In 
punched  into  the 


reduction  system  using  the  Hewlett-Packard  Model  20 
calculator  can  be  easily  adapted  to  fully  automated 
this  case,  the  raw  data,  P,  would  not  be  mcuiually 
calculator,  but  it  would  be  directly  coupled  in  from 
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AURORA  Fncility 


DOSIMETRY  DATA  ^or  Shot  1757 


100  ns/in. 


Date:  4 Mar  75 

Zero  time:  1100 

Area:  194  in. -ns 

PeeOc  Height:  1.59  in. 
Eff.  V/idth:  122  ns 


Figure  17 . 


Dosimeter  No. 

DOSE 

(rads (Si)) 

DOSE  RATE  peak 

(rads(Si)/s) 

1 1 

1.632E  04 

1.339E  11 

12 

1 . 507E  04 

1.236S  11 

13 

1 . 705E  04 

1.399E  11 

14 

1.600E  04 

1.312E  11 

15 

1.659'^  04 

1.361E  11 

16 

I.5IAE  04 

1.24^  ' 11 

17 

1.597E  04 

1 .3IOE  11 

18 

1 .661E  04 

1 .362':  1 1 

Example  of  dose  and  dose-rate  data  sheet. 


the  interface  unit  by  an  appropriate  data  coupler.  Alternatively,  the 
data  could  be  recorded  on  a magnetic  tape  cassette  and  later  played  back 
into  the  Model  20.  Such  an  automatic  data  recording  system  using  an 
Electronic  Processors  STR-200I  incremental  cassette  recorder,  a 
parallel“tO“Serial  converter,  and  the  HP  11202A  TTL  I/O  Interface  Card 
is  now  under  construction. 

Combining  such  an  automatic  data  reduction  system  with  an 
automatic  TLD  changer  could  result  in  complete  hands-off  reading  and 
data  reduction  of  up  to  200  TLD*s  at  a time. 


7 . CALIBRATION 

Two  methods  are  routinely  used  to  calibrate  and  monitor  the 
calibration  of  the  AURORA  TLD  readers. 

(a)  A number  of  TLD's  are  exposed  for  a fixed  time  interval  in  a 
Co  source  whose  exposure  rate  is  well  known.  The  corresponding  dose 
deposited  by  AURORA  radiation  can  be  calculated  and  compared  with  the 
dose  calculated  from  the  TLD  reader  output. 
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(b)  A number  of  TLD's  is  exposed  to  an  AURORA  pulse  alongside  an 
absolute  silicon-dose  calorimeter.  The  silicon  dose  as  calculated  from 
the  reader  output  is  compared  with  the  dose  measured  by  the  silicon 
calorimeter. 


7.1  Cobalt-60  Calibration 

The  success  of  Co^^  calibration* *  depends  on  two  things.  First, 
the  exposure  rate  at  a given  position  in  a Co^O  source  must  be 
accurately  known.  The  rate  is  usually  determined  with  air-equivalent 
ionization  chambers  calibrated  by  the  National  Bureau  of  Standards.  Dhe 
second  is  the  ability  to  convert  Co^®  exposure  in  roentgens  to  AURORA 
silicon  dose  in  rads (Si).  This  ability  requires  a knowledge  of  the 
AURORA  spectrum,  which  we  do  not  have.  However,  some  careful  estimates 
yield  results  that  are  believed  to  be  quite  good. 


The  assumption  behind  this  derivation  is  that  the 
thermolximinescence  response  of  a TLD  will  be  the  same  for  equal  dose 
(rads  (CaF2  :Mn) ) , whether  this  dose  is  deposited  by  Co^^  or  by  AURORA 
bremsstrahlung . 


First,  one  calculates  „ , deposited  by  Co^®.  To 

CaF2 

one  exposes  the  TLD  in  a G.43-cm-thick  aluminum  capsule, 
unperturbed  free  field  exposure,  Xq , in  the  Co^^  source, 
exposure,  X,  at  the  TLD  is 


enstire  CPE 
Given  the 
the  true 


X = Xoexp  = 0.971X0 

where  is  the  mass  energy  absorption  coefficient  for  Co^O  radiation 

in  aluminum,  and  t is  the  aluminum  thickness. 


To  determine  D from  X,  one  assxames  that  the  TLD  in  the 
CaF2 

aluminum  capsule  behaves  like  a Bragg-Gray  cavity.^  That  is,  primary 
photon  interactions  occur  only  in  the  aluminum,  and  all  of  the  dose 
deposited  in  the  CaF^  is  by  secondary  electrons  from  the  al\iminum.^ 


^F.  H.  Attix,  W.  C.  Roeschf  and  F. 
2nd  Ed.,  Academic  Press,  New  York  (1966). 

*F.  H.  Attix  provided  this  analysis. 

^This  assumption  is  not  necessarily 
thickness  of  the  TLD;  nevertheless,  it 
stopping  power  ratio  approximately  equals 
coefficient  ratios. 


Tochilin,  Radiation  Dosimetry , 


valid,  because  of  the  finite 
is  justified  here,  because  the 
the  mass  energy  absorption 


(26) 
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From  the  definitions  of  the  rad  and  the  roentgen,  we  know  that, 
for  Co^^  radiation, 


D . = 0.869X  . 

air 


(27) 


Considering  just  the  primary  photon  interactions  in  aluminum. 


= D . = 0.966D  . 


(28) 


(values  of  y /p  are  found  in  standard  tables^' If  all  the  dose 
en 

in  the  CaF2  comes  from  secondary  electrons  produced  in  the  aluminum, 
the  average  energy  of  these  secondary  electrons  would  be  approximately 
300  keV.  Then, 


^CaF2 


sCaF2 
300  key 
A1 

^300  keV 


D 


A1 


1.013D 

A1 


(29) 


(values  of  S are  found  in  standard  tables^' where  is  the  electron 

E 

stopping  power  of  material  Z for  electrons  of  energy  E.  Combining 
equations  (27)  to  (29)  gives 


D = 1.013  X 0.966  X 0.869X  = 0.850X  (30) 

CaF2 

for  the  dose  deposited  in  CaF^  by  Co^^  gamma  radiation. 


^F.  H.  Attix , W.  C.  Roesch,  and  E,  Tochilin , Radiation  Dosimetry , 
2nd  Ed.,  Academic  Press,  New  York  (1966). 

^Ann  T.  Nelms,  Energy  Loss  and  Range  of  Electrons  and  Positions , 
National  Bureau  of  Standards  (U.S.),  Circ.  577  (1956). 

Hubbell,  Photon  Cross  Sections,  Attenuation  Coefficients , and 
Energy  Absorption  Coefficients  from  10  keV  to  100  GeV , National  Bureau 
of  Standards  NSRDS-NBS29  (August  1969). 

Berger  and  S.  Seltzer,  Tables  of  Energy  Losses  and  Ranges  of 
Electrons  and  Positrons , National  Aeronautics  and  Space  Administration 
NASA  SP-3012  (1964) . 
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At  AURORA,  photon  attenuation  experiments  indicate  that  the 

average  photon  energy  is  between  1 and  2 MeV.  One  wishes  to  find  the 

dose  in  CaF2,  given  a certain  silicon  equilibrium  dose,  D ..  Also,  one 

Sx 

ensures  that  one  has  CPE  by  enclosing  the  TLD  in  a sufficiently  thick 
aluminum  capsule.  Assuming  CPE, 


D 


A1 


/ / \ 

TMeV  „ 

/ Si 

MeV 


0.966D^.  . 
Si 


(31) 


Again  assuming 
cavity. 


that  the  aluminum-CaF2  system  approximates  a Bragg-Gray 


CaF2 


gCaFa 
600  keV 
A1  °A1 

600  keV 


1.012D,,  . 
A1 


(32) 


Combining  equations  (31)  and  (32), 


D = 1.012  X 0.966D  . = 0.978D^.  (33) 

CaF2  Si  Si 

for  the  dose  deposited  in  CaF2  by  AURORA  bremsstrahlung. 

Now,  with  the  assumption  that  the  TL  response  of  the  same  TLD 
is  equal  for  equal  absorbed  dose  regardless  of  the  source,  exposure  to 
0.850  R(Co^O)  will  give  the  same  light  output  as  the  absorption  of 
0.978  rads (Si) (AURORA) ; or 

0.850  X(Co®°)  = 0 . 978D  (AURORA)  (34) 

1 


whence 


D .(AURORA)  = 0.870X(Co^°) . 
Si 


(35) 


This  equation  (35)  is  our  basic  calibration  relationship. 
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By  use  of  the  calibration  relationship,  a whole  calibration 
curve  like  that  in  figure  16  can  be  established.  Once  it  is,  however, 
it  need  not  be  repeated,  because  its  shape  does  not  change,  even  from 
one  TLD  production  lot  to  the  next.  To  ensure  that  the  reader  is 
maintained  in  calibration,  it  is  necessary  merely  to  recheck  one  point 
on  figure  16  periodically  as  follows:  A representative  sampling  of 
aurora's  current  lot  of  TLD's  is  exposed  in  a Co^^  source  to  an 
unperturbed  exposure  of  2.0  x 10^  R.  The  true  exposure  at  the  TLD  is 

0.971  X 2 X 10^  = 1.94  x lo^  R (eq.  (26)). 


This  exposure  corresponds  to 

0.870  X 1.942  X lo^  = 1.69  x lo^  rads(Si)  (eq.  (35))  . 

Reference  to  figure  16  reveals  that  for  P = 1.70  x lo^ , D/P  = 0.994,  and 
D = 1.69  X 10^  rads(Si).  So  TLD's  exposed  to  an  unperturbed  exposure  of 
2 X 10^  R should  yield  a reading  of  exactly  1.70  x at  1 hr  after 

exposure. 


The  calibration  of  the  TLD  readers  is  checked  periodically, 
sometimes  every  day,  by  reading  10  of  these  exposed  TLD's,  correcting 
for  fading,  and  comparing  the  mean  of  the  10  readings  to  1.70  x lO^. 
Since  10  readings  are  a relatively  small  sample  of  the  population,  the 
mean  P determined  this  way  is  typically  different  from  1.70  x 10^.  To 
test  for  the  statistical  significance  of  this  variance,  the  Student's 
t-test  is  used.^^  This  is  not  nearly  as  laborious  as  it  sounds,  since 
the  t-test  for  the  significance  of  a sample  mean  is  calculated 
automatically  with  the  HP  Model  20  calculator.  If  the  t-test  reveals  a 
significant  difference  between  P and  1.70  x 10^,  the  sensitivity  of  the 
TLD  reader  is  adjusted  to  bring  it  back  into  calibration. 

7.2  Silicon-Calorimeter  Calibration 


Calibration  by  the  silicon  calorimeter  dispenses  with  all  the 
assumptions  made  in  the  previous  paragraphs  and  compares  the  TLD  reader 
response  directly  to  an  absolute  silicon  dose  standard,  irradiated  in 
the  AURORA  environment.  *' Absolute"  means  that  the  silicon  calorimeter 
need  not  be  calibrated  against  a known  radiation  standard.  It  is 
necessary  merely  to  be  able  to  measure  the  change  in  the  resistance  of  a 
thermistor  very  accurately  to  determine  the  silicon  dose.  The 
construction  and  operation  of  this  calorimeter  is  fully  described 


^^William  Volk,  Applied  Statistics  for  Engineers,  2nd  Ed.,  McGraw-Hill 
Book  Co.,  New  York  (1958). 
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elsewhere.^  Briefly,  it  consists  of  a tiny  thermistor  embedded  in  a 
silicon  pellet  that  is  surrounded  by  a silicon  mantle  thick  enough  to 
guarantee  CPE.  The  whole  calorimeter  is  thermally  isolated  from  its 
surroundings. 

The  rise  in  temperature  of  the  silicon  pellet  upon  irradiation 
depends  only  on  the  specific  heat  of  silicon  and  the  absorbed  dose.  The 
temperature  rise  is  determined  by  the  thermistor,  whose  (negative) 
temperature  coefficient  of  resistance  is  well  known.  Measurements  of 
Rq,  the  equilibrium  resistance  of  the  thermistor,  and  AR,  the  change  in 
resistance  upon  irradiation,  are  the  only  measurements  required  to 
calculate  the  dose.  The  resistance  parameters  are  measured  with  a good 
Wheatstone  bridge  (Keithley  Model  520) . The  dose  per  unit  resistance 
change,  D /Ar,  is  a nonlinear  function  of  Rq.  To  facilitate  this 
calculation,  this  function  was  approximated  by  a fifth-degree 
polynomial,  and  the  dose  calculation  is  done  on  the  calculator. 


For  the  typical  dosimeter  calibration,  eight  TLD's  are  placed 
close  to  the  calorimeter  (fig.  18) , and  the  whole  assembly  is  exposed  to 
an  AURORA  pulse.  Measurements  for  AR,  Rq,  and  Pj  to  P9  for  the 
dosimeters  are  put  into  the  calculator.  The  calculator  calculates  first 
the  dose  as  measured  by  the  calorimeter  and  then  the  dose  corresponding 
to  Pi  to  Pa  using  equation  (22).  It  then  does  a t-test  to  test  whether 
any  variance  between  D (calorimeter)  and  5 (TLD)  is  statistically 

O X wX 

significant. 

When  a TLD  reader  is  calibrated  by  the  Co^^  method,  the  silicon 
calorimeter  test  indicates  also  that  the  reader  is  calibrated.  In  other 
words,  the  two  methods  do  not  differ  significantly.  Both  methods  are 
therefore  equally  valid  for  calibrating  the  TLD  readers,  and  both  are 
used  routinely. 


^F.  H.  Attix,  Dosimetry  for  the  Aurora  Facility , Final  Report,  Naval 
Research  Laboratory , Washington,  DC,  NRL  Memorandum  Report  2552,  Revised 
(March  1973). 
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(negative  No.  099-74) 

Figure  18.  Absolute  silicon  calorimeter. 

8.  PRECISION  AND  ACCURACY 
8.1  Definitions 


When  a given  measurement  is  repeated  a number  of  times,  the 
resulting  values  will,  in  general,  not  agree  exactly.  These 
discrepancies  are  caused  by  experimental  errors.  A measurement  with 
small  experimental  errors  has  high  precision. 

If  all  measured  values  are  in  error  by  the  same  amount,  this  is 
a systematic  error.  Measurements  with  small  systematic  errors  have  high 
accuracy. ^ ^ 

A good  measure  of  the  experimental  error  is  the  standard 
deviation.  If  the  experimental  errors  are  random,  a large  number  of 
measurements  of  the  same  quantity  will  have  a normal  distribution.  The 
standard  deviation  of  a number  of  normally  distributed  measurements,  xi, 
can  be  estimated  by 


^^Yardley  Beers,  Introduction  to  the  Theory  of  Error,  Addison-Wesley 
Publishing  Co.,  Inc.,  Reading,  MA  (1956). 
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(36) 


^ S x2  - Nx2jis 

where  x is  the  arithmetic  mean  and  N is  the  number  of  measurements.  The 
probability  of  a measurement  differing  from  the  mean  by  more  than  one 
standard  deviation  (1  is  0.32;  the  probability  of  a 2-a  variance  is 
0.045. 


8.2  Precision  of  TLD  Measurements 


The  percentage  standard  deviation  of  a large  number  of  TLD 
measurements  of  the  same  dose  has  been  experimentally  determined;^ 


TLD  %s(x) 


Teflon  disc  3.15 
Solid  CaF2:Mn  5.4 

This  is  a measure  of  the  precision  of  dose  measurements  made  with  these 
TLD's.  What  this  means  is  that  95  percent  of  all  dose  measurements  are 
within  ^6.3  percent  of  the  mean  for  Teflon  disc  measurements  and  within 
±11  percent  for  solid  CaF2:Mn  measurements. 

8.3  Accuracy  of  TLD  Measurements 

The  cimount  by  which  the  mean  of  a large  number  of  measurements 
differs  from  the  "true"  value  of  the  quantity  being  measured  is  a 
measure  of  the  accuracy  of  the  measurement.  The  accuracy  of  TLD 
measurements  of  dose  is  difficult  to  estimate  since  one  does  not  know,  a 
priori,  the  systematic  errors  in  the  measurements.  Two  independent 
measurements  of  the  dose  (i.e. , use  of  TLD's  calibrated  by 
Co^^  exposure  and  use  of  the  silicon  calorimeter)  generally  agree  to 
better  than  3 percent.  From  that  fact,  one  can  guess  that  the  accuracy 
of  these  measurements  is  better  than  5 percent. 
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